Laser-induced currents along molecular wire junctions 
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The treatment of the previous paper is extended to molecular wires. Specifically, the effect of 
electron-vibrational interactions on the electronic transport induced by femtosecond ui + 2u laser 
fields along unbiased molecular nanojunctions is investigated. For this, the photoinduced vibronic 
dynamics of irans-polyacetylene oligomers coupled to macroscopic metallic leads is followed in a 
mean-field mixed quantum-classical approximation. A reduced description of the dynamics is ob- 
tained by introducing projective lead-molecule couplings and deriving an effective Schrodinger equa- 
tion satisfied by the orbitals in the molecular region. Two possible rectification mechanisms are 
identified and investigated. The first one relies on near-resonance photon-absorption and is shown 
to be fragile to the ultrafast electronic decoherence processes introduced by the wire's vibrations. 
The second one employs the dynamic Stark effect and is demonstrated to be highly efficient and 
robust to electron-vibrational interactions. 



I. INTRODUCTION 

Considerable effort has been devoted to studies of the 
properties of molecular wires [l|, 0, d, 0, H, H 0| , since 
they constitute natural candidates for the ongoing minia- 
turization of electronic devices. Generally, the focus is 
placed on the transport properties of metal-molecule- 
metal junctions subject to a bias voltage. In this regime 
the metallic leads are the main source of electrons for the 
transport while the molecular system serves as a trans- 
porting medium that can be chemically functionalized to 
modify the I-V characteristics of the junction. 

Here we consider an alternative situation in which the 
junction is not subject to a bias voltage and where the 
molecule serves as the main source of transporting elec- 
trons. The composite system is taken to be spatially 
symmetric and net currents are induced using lasers with 
frequency components u and 2uj. Such fields give rise to 
phase-controllable transport in symmetric systems even 
when they have a zero temporal mean [8|, |9| . This rec- 
tification effect first appears in the third order response 
of the system to u> + 2oj fields. At this order the sys- 
tem mixes the frequencies and harmonics of the incident 
radiation generating a phase-controllable zero-frequency 
(DC) component in the response. The setup is of interest 
since it can produce ultrafast currents and may lead to 
the development of molecular switches that operate on a 
femtosecond timescale. 

When using uj + 2uj fields to induce electronic trans- 
port it is crucial to consider the influence of the sys- 
tem's vibrations on the laser scenario. This aspect 
becomes particularly relevant in molecular wires since 
molecules nearly always exhibit vibronic couplings that 
can substantially modify the photophysics of the junc- 
tion and the effectiveness of the scenario. This charac- 
teristic of molecular electronics marks an important dif- 



ference relative to the behavior of rigid semiconducting 
solids, and is currentl y th e focus of vigorous investigation 
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In this paper we characterize the influence of vibronic 
couplings on rectification induced by femtosecond uj + 2uj 
laser pulses along unbiased molecular wires. As a sys- 
tem we consider a irans-polyacetylene (PA) oligomer 
connected by its ends to macroscopic metallic leads. 
The leads are treated as rigid semi-infinite tight binding 
chains and the oligomer is described via the Su-Schrieffer- 
Heeger (SSH) Hamiltonian [U,[ig[. 

This work is a continuation of our studies in Ref. [20 in 
which we considered the problem of laser-inducing rec- 
tification along isolated PA chains. Naturally, it shares 
many of the challenges encountered therein. In partic- 
ular, the problem of how to effectively exert laser con- 
trol on the electronic dynamics in the presence of lattice- 
induced decoherence persists. However, it differs signif- 
icantly from the case of isolated chains since the cou- 
pling to the metallic contacts introduces a decay route 
for energetic electrons and adds an additional element of 
complexity to the theoretical description. 

Besides offering insight into vibrational effects in 
molecular electronics, this study contributes to the 
emerging direction of using lasers to modify the trans- 
port properties of nanojunctions in which several promis- 
ing applications are being recognized. For instance, it 
has been demonstrated that light can induce a switch- 
ing behavior in the conduction when, upon photoex- 
citation, the molecule undergoes substantial conforma- 
tional change [HI [22l H3- ^ n addition, recent theoret- 
ical work suggests that lasers may enhance [HI, HE] or 
suppress [H, [27| transport induced by a bias or, alter- 
natively, gene rate currents in unbiased but asymmetric 
junctions [23, El HI- 

Additional motivation for this study was provided by 
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a rather puzzling observation made in Ref. [30(. In that 
work, the authors analyzed vibrational effects on laser- 
induced rectification along nanojunctions. The model 
presented there refers to a two-site molecular wire driven 
by a continuous wave ui + 2w field. In it, each site is 
coupled perturbatively and independently to one local 
bath of harmonic oscillators in thermal equilibrium, and 
an Ohmic spectral density is adopted to describe the 
phonon-induced relaxation. This work concludes that 
vibronic couplings enhance the current induced by the 
symmetry breaking field by more than one order of mag- 
nitude. This observation seems to be in contradiction 
with our results in Ref. [2^ in which the vibronic couplings 
were shown to have strong detrimental effects on the con- 
trol, specially when long pulses are employed. Further, 
no insight into this rather nonintuitive behavior was pro- 
vided. 

In addition to clarifying this observation, and faith- 
fully characterizing vibrational effects on laser rectifica- 
tion along molecular nanojunctions, a novel mechanism 
to induce ultrafast currents along molecular wires that re- 
lies on the dynamic Stark effect instead of near-resonance 
photon absorption is discussed here [3l[. As shown be- 
low, this mechanism is remarkably robust to electron- 
vibrational couplings and is able to induce large currents 
in molecular wires with efficiencies > 90%. 

As in Ref. [20l . the photoinduced vibronic dynamics 
of the wire is followed in a mean-field mixed quantum- 
classical approximation \3^, [H, [U, H|| . Electronic deco- 
herence and relaxation is incorporated by averaging over 
an ensemble of quantum-classical trajectories with initial 
conditions obtained by importance sampling the ground- 
state nuclear Wigner phase space distribution. In this 
manner, the simulations take into account the fact that 
the vibrational modes in typical molecular wires are non- 
local and constantly kept out of thermal equilibrium ei- 
ther through direct interaction with the field or through 
exchange of energy with the driven electrons. 

A reduced description of the dynamics of the composite 
system is obtained by deriving an effective Schrodinger 
equation satisfied by the orbitals in the molecular region. 
As shown, in the wide bandwidth limit the lead regions 
are mapped into negative imaginary (absorbing) poten- 
tials. In addition, we employ projective lead-molecule 
couplings that incorporate basic aspects of the Fermi sea 
while avoiding the necessity to explicitly follow the dy- 
namics of the essentially infinite number of electrons in 
the leads. This model is only appropriate to describe the 
short time dynamics of unbiased junctions since, in it, 
charge is not permitted to flow from the leads into the 
oligomer. 

This manuscript is organized as follows: In Sec. ITT1 we 
briefly describe the model (Sec. Ill A|) . and determine the 
effective mean-field equations of motion for the nanojunc- 
tion (Sec. Ill Bj) . Subsequently, we derive an expression for 
the photoinduced currents ( Sec. Ill Cj) that solely depends 
on molecular properties. Our main results are stated in 
Sec. Mil Specifically, in Sec. IIII Al we study the scenario 



in its usual regime where the laser frequencies are tuned 
at or near resonance with one of the electronic transition 
frequencies. Subsequently, in Sec. IIII Bl we introduce the 
Stark shift rectification mechanism. The emerging trends 
are summarized in Sec. IIVI 



II. MODEL AND METHODS 
A. Model 

As a system we consider a irens-polyacetylene oligomer 
connected by its ends to macroscopic metallic leads. 
Since typical metals screen electric fields that have a fre- 
quency below the plasma frequency, so that electromag- 
netic radiation from the optical or the infrared spectral 
range is almost perfectly reflected at the surface [3y|, we 
assume that the radiation field only interacts with the 
oligomer chain. 

The composite system is described as a one- 
dimensional lattice in which each site n corresponds to 
the position of an atom, and is defined by the Hamilto- 
nian 

H(t)=H L + H s _ L + H s (t) + H s _ R + H R . (I) 

The left (L) and right (R) leads are treated as semi- 
infinite rigid tight-binding chains with hopping param- 
eter tlead, SO that 

■Hl = — ^lcad y ] y ] ( C n+l.s C n,s + c n.s c n+l,s), 

n<0s=±1 , (2) 

Hr = -<lcad ^2 ( C l+l,s C n,s + cl l s C n+hs ). 

n>Ns=±l 

Here cjj s (cn, s ) creates (annihilates) a fermion in site n 
with spin s and satisfies the usual fermionic anticom- 
mutation relations {c n ,ST^ msl } = 5 n , m 8 SyS i . The N- 
membered oligomer chain with Hamiltonian H^{t) — 
H^(t) + i?g h (i) is situated between sites n = 1,...,N 
and is described by the SSH model [Ig, [T9| coupled to a 
laser field E(t) in dipolc approximation. Specifically, the 
electronic part of Hg is 

N-l 

H s(t) = Mo + a(u n+1 - w„)](<4 +M c„, s 

n—l,s , . 

3 

N y ' 

~^ c n,s c n+l,s) 4" \e\ ^ ^ X n cl l s C ntS E(t), 
n=l,s 

where io is the hopping integral for zero displacement, 
u n is the monomer displacement of site n from its per- 
fectly periodic position na, and a describes the electron- 
ion coupling between neighboring sites. In turn, x n = 
(na + u n ) is the position operator for site n, a the lattice 
constant, and — |e| the unit electronic charge. The lattice 
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is described by 



N 



Pn , K 



N-l 



N 



2M 2 

n—1 n— 1 



£ (« n +i - u n ) 2 ~\e\ £ x n E(t), 



(4) 

with force constant i-C and (CH) group mass M, and 
where p n is the momentum conjugate to u n . Last, the 
coupling between the oligomer and the leads 



H S -l = -tcoup £ ( c i, s c o, s + 4, s c m) 

s=±l 

Hs-R = —tcoup £ (c]v +1 s Cjv,s + Cjy^CiV+l^), 



(5) 



s=±l 



is taken to be at a single site only with coupling constant 

tnonn • 



B. Equations of motion 

The electron-vibrational dynamics of the composite 
system in the presence of a radiation field is followed in 
the mean-field (Ehrenfest) mixed quantum-classical ap- 
proximation introduced in Sec. IIB of Ref. 20. As before, 
the nuclei satisfy classical trajectories governed by 



u n (t) 

Pn(t) 



_ Pn{t) _ 
M ' 

= - K (2u n (t) - u n+1 (t) - u„_i(t)) 

Pn,n— 1 

(*)} 

-\e\E(t) (p n , n (t)-l), 



(6) 



where p n ,m(t) = 2 s (*(*)l c n, s c m>*l*(*)) is the reduced 
electronic density matrix. The chain is taken to be 
clamped so that u±(t) =upf(t) =0 and pi(t) = pjv(i) = 
0, and Eq. © is valid for n = 2, • • • , N - 1. In 
turn, the many electron wavefunction ^(t)) satisfies the 
Schrodinger equation and, hence, the dynamics of p n>m (t) 
is given by 

ih^- Pn , m (t) = 5^(*(t) I [4,.Cm,», H e i ec (t)} (7) 



where 



S-R 



(8) 



is the electronic part of the Hamiltonian. For an initial 
electronic state |^(0)) that is well described by a single 
Slater determinant, the reduced density matrix admits 
the orbital decomposition 

Pn,m 

W = H^ 6 ^)' s l n ' s >( m ' s\e(t),s)f(e, s) (9) 



where /(e, s) is the time-independent initial distribution 
function that takes values or 1 depending on the oc- 
cupation of each level with energy e and spin s, and 



|n, s) — cj, JO) where |0) represents the vacuum state. In 
order for p n<m (t) to satisfy Eq. 0, the orbitals |e(t),s) 
must be solutions to the single-particle time-dependent 
Schrodinger equation 



ih— \e(t),s)=H elec (t)\e{t),s) 



(10) 



with initial condition \e(t = 0),s) = \e,s), where |e, s) 
denote the eigenorbitals of -ff c icc at preparation time. 

In principle, by integrating the above equations self- 
consistently it is possible to follow the coupled dynam- 
ics of both electronic and nuclear degrees of freedom in 
the composite system. In practice, the problem is in- 
tractable since the number of equations to be integrated 
is unbounded due to the large number of electrons in 
the leads and the spatial unboundedness of the compos- 
ite system. Further, the simplifications that have been 
developed for periodic systems are not applicable here. 
A reduced description of the dynamics is obtained by 
first deriving, in Sec. Ill B li an effective single-particle 
Schrodinger equation that is satisfied in the molecular re- 
gion. This eliminates the necessity of integrating Eq. (flU]) 
over the complete spatial domain of the composite sys- 
tem. Subsequently, in Sec. Ill B 21 we specialize the cou- 
pling between the molecule and the leads using projec- 
tion operators so that basic aspects of the lead-molecule 
coupling are captured without explicitly following the dy- 
namics of the electrons originally in the leads. Then, in 
Sec. Ill B 31 we derive an expression for the memory ker- 
nel of the resulting equations, which is then employed in 
Sec. lIIB4l to obtain the Markovian limit of the dynamics. 



1. Effective non- Markovian Schrodinger equation 

Henceforth, we drop spin labels since orbitals with the 
same energy but opposite spin satisfy the same equation 
of motion under H e i ec (t). We begin by partitioning the 
orbitals that form p n ,m(t) into the three spatial regions 
of the composite system 



\e(t)) = |e L (i)> + |e s (t)> + \e R (t)). 



(11) 



Here \ep(t)) — X)nep( n l e W)l n ) describes the part of the 
orbital in region P, where P can be the left lead (L), 
right lead (R) or system (S). In terms of this partition 
the single-particle Schrodinger equation (fTOf reads 



ih— 
dt 



>L(*)> 




ks(i)> 




|e R (*)> 













<*(*)> 
es(*)) 

(12) 

where we have explicitly allowed for time-dependence in 
the lead-molecule coupling since this characteristic will 
arise when introducing projection operators into the dy- 
namics. 

We seek an effective equation of motion for | eg) in 
which it is not necessary to explicitly follow the dynamics 
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of |£l) or |ejt) [33|- This can be accomplished by Laplace 
transforming the evolution equation for the lead compo- 
nents, solving the resulting algebraic equations and in- 
verse Laplace transforming by taking into account the 
convolution theorem. The procedure yields: 

Mt)) = U?{t)\ep{Q)) + \ f U^t-r)H s _ p {T)\e s {r))AT, 

(13) 

where (3 — L or R denotes the left or right lead, respec- 
tively. The first term in Eq. (|13|) describes the evolution 
of the lead components of the orbitals in the absence 
of coupling to the molecular system. Naturally, the dy- 
namics is determined by the evolution operator for the 
isolated leads U@(t) = exp(—iH^t/h). In turn, the in- 
fluence of the lead-molecule coupling is mapped into a 
convolution integral with U°(t — r) acting as the mem- 
ory kernel. 

By introducing Eq. (| 13|) into the evolution equation for 
|es) we arrive at an effective non-Markovian Schrodinger 
equation for the system part of the orbital: 

ih±\e s (t))=H 3 (t)\e s (t)) + [#s-^)^WM0)> 

/3=L, R. 

+ / Uf(t - r) J ff s ^(r)|e s (r))dr . 

in j j 

(14) 

The first term corresponds to the single-particle 
Schrodinger equation for the isolated system. The re- 
maining terms describe the transfer of population ampli- 
tude between the system and the leads. 

The advantage of Eq. (HU) over Eq. (jlOp is that it per- 
mits following the dynamics of the composite system by 
integrating the equations of motion in a finite spatial re- 
gion. The influence of the leads has been mapped into 
a term that is nonlocal in time and a residual term de- 
pendent on the initial state of the leads. Equation (fTi]) 
is also a good starting point for various approximations. 
In particular, the Markovian limit arises naturally when 
the energy bandwidth of the metallic leads is much larger 
than any other characteristic energy of the system. 

2. Projective molecule-lead coupling 

Solving the above integro-diffcrential equation 
[Eq. (fT4"|) ] exactly for all the initially occupied levels in 
the composite system is a problem as formidable as the 
original one, and an effective way to introduce the leads 
into the dynamics is required. Physically, the contacts 
block any electrons from the molecule with energy less 
than the Fermi energy and act as a thermal source of 
electrons to the molecular system. We are interested in 
the case where there is no bias voltage across the bridge 
and in short simulation times that limit the degree of 
relaxation of electrons from the leads to the molecular 
states. In this regime, the main effect of the leads is to 



accept electrons that have energy higher than the Fermi 
energy ep and to block them otherwise. We model this 
property using projection operators that restrict the 
interaction between the molecule and the leads to only 
those molecular levels that are above the Fermi energy. 
Consider a projection operator (V 2 = V) 

V(t) = V s (t) +V L +V R (15) 

where 

^=Emm, as) 

with j3 —L or R, projects onto single-particle states of 
the left/right lead. In turn, 

r s (t)= E l7sW)(7s(t)| (17) 

is a projection operator onto the instantaneous light- 
dressed eigenorbitals |7s(*)) of the molecular system that 
reside above the Fermi energy. The latter are defined by 
the eigenvalue relation 

tfs ol (i)l7sW)=6 7 W| 7 s(i)), (18) 

where HI ! (i) is the electronic Hamiltonian of the oligomer 
at time t given in Eq. ©. The time-dependence of |7s(i)) 
[and V s (t)] is a result of the time-dependence in the elec- 
tronic Hamiltonian introduced by the field and the nu- 
clear motion. An evident but important property of the 
projection operator is that V^V® — V^V 5 = 0. 

We now wish to use V to restrict the lead-molecule in- 
teraction. The projection operator cannot be introduced 
directly into the second quantized version of the lead- 
molecule coupling [Eq. (J5J], Rather, it is necessary to 
modify the first-quantized interaction Hamiltonian and 
then second quantize it. In first quantization, the origi- 
nal coupling takes the form -ffs-,3 = Y^iLi n s-p(i), where 
the label % refers to the coordinates of the i-th electron 
and the sum runs over all M electrons in the system. 
The nearest- neighbor couplings of Eq. ([5]), give rise to 
the following matrix elements of /is-/3, 

(n\h s . L \m) = -i C ou P (^n4^m,o + <5n,o<W), 
{n-\hs-R\m) = —t conv (5 n ^5 m .N+i + S n ,N+i5 m ,N)- 

By introducing V into the first-quantized coupling Hamil- 
tonian, so that 

M 

only those electrons in the desired subset of electronic lev- 
els are allowed to move from the molecule to the leads, 
and vice versa. The modified coupling operator is Hermi- 
tian, and thus diagonalizable, and hence can be cast into 
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second quantized form using the usual procedure [38| , to 
obtain 

Hs-h(t) = -icoup ^2 Vn,o(t)Pl,rn(t)cl s C m ,s + H.C, 
Hs-n(t) = — t coup ^ 'P n ,iv(^)^ :) 7V+l,m(^)c^ s C m!S + H.C, 

(21) 

where 7-Vi, m = (n\V\m) and H.c. stands for Hermitian 
conjugate. 

We specialize our considerations to the case in which 
electrons in the molecule are allowed to leak into the leads 
but where no electrons that are originally in the leads can 
enter the molecule. This is modeled by treating the leads 
as electron-less and by coupling the excited electrons in 
the molecular system to every level in the leads. Since 
we now project onto a complete lead basis V n fi = #n,o 
and Vn.N+i = S n ,N+i, and the lead-molecule interaction 
terms assume the simplified form: 

H S . L (t) = -t coup Y, E KnitH^s +H.C, 

nes s=±i 



nSS s=±l 



(22) 



The leads and molecule are taken to be initially de- 
tached and, since the projective coupling Eq. (|2"2"j) al- 
ready incorporates the Fermi blockade imposed by the 
leads, we only consider orbitals that are initially in the 
molecular region so that \ep(t = 0)) = for the states of 
interest. Taking into account the above, and introducing 
the projective coupling ([22]) into Eq. (fl4|) . it follows that 

i N 
ih-(n\e s (t)) = ]>>|tf s cl (i)|m)(m| £s (t)) 



m—1 

,2 N rt 



m— 1 

+ -PnAt)U% + i,N + i(t ~ r)V s N . m (r)} (m|e s (r))dr. (23) 



3. The memory kernel 

To proceed further an explicit expression for the mem- 
ory kernels in Eq. (|23|) is required. The task is to find 
the matrix elements of the evolution operator U@(t) — 
exp(—iHpt/K) for the left (f3 =L) and right (f3 =R) semi- 
infinite tight binding chains in site representation. For 
this, first note that the lead Hamiltonian [Eq. ([2])] can be 
diagonalized by performing the basis transformation: 



o \ 1/2 n 

h) S 

7 fe=i 



sin [(n - n^kir/tt] c fc , s , (24) 



where f2 is the number of sites in the leads and np is the 
molecular site that is connected to lead /3, i.e. = 1 
and ur = N. Inserting Eq. (JUJ into Eq. ^ yields 



Hp = -2<i oad ^ Y cos(fc7r/n)4 



fc,s C *V 



(25) 



k=l s=±l 



where we have exploited the fact that for large f2, S^y — 

h E«=i ex P M fc - k')ir/n). 

The matrix elements of the evolution operator in site 
representation n are then given by: 



U£ m (t) = (n\uP(t)\m) 



^2 s i R [i n ~ np)kir/Q] 



k=l 



sin [(to — np)kir/£l\ exp [2it\ ca< \ cos (kn/Q)t/h], 

where we have used Eqs. (|24l) and (|23|) . For large SI the 
sum can be approximated by an integral to give: 



1 

71" ./,. 



{ cos [(n + m — 2np) 



cos 



[{n — m)0] } exp [2rfi ea( j cos(8)t / 'K\ AO 



_ -n-rn j ( -^lcad 

■ I J n —m, I . t 



V h 



fJ n+m—2n(3 \ t- b I i 



(26) 



where J„(z) = '^-^ f£ e zzcos9 cos(n9) AO is a Bessel func- 
tion of the first kind of order n [39| . We note that U% m 
satisfies the equations of motion, as well as the boundary 
and initial conditions: 



h-Ul m (t) =J2(n\Hp\r)U r e m (ty, 



rG/3 



U^, m (t) = U^ np (t) = 0; tf£ m (t = 0) = S n , m . 
Of particular relevance is 

2Ji (*2%*Mf) 

K{t) = uUt) = U% +1 , N+1 (t) = \\J t (27) 



where we have used the identities J_„(z) = (— l) n J n (z) 
and J„_i(z) + J n +i(z) — — J n (z). This quantity deter- 
mines the kernel of the convolution integral in Eq. (|23|) . 
It is characterized by a pronounced peak around t = 
followed by oscillations that decay asymptotically as 
i / n \ 3/2 

^5= ( fi dt J ■ This asymptotic dependence is exploited 
in the next section to derive the Markovian limit of the 
equations of motion. 
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4- Wide bandwidth approximation 

We now invoke the wide bandwidth approximation in 
which the energy bandwidth of the leads (4ti oa( j) becomes 
the dominant energy in the problem and, consequently, 
h/tiead the fastest time scale. In this regime, the memory 
kernel ]C(t — r) becomes sharply peaked around r = t and 
the orbitals do not change appreciably on the time scale 
in which JC(t — r) varies. It follows that the convolution 
integrals in Eq. (|23[) can be approximated by 



by Eq. ([6]). The orbitals determining p n , m (t) [Eq. {9j] are 
the initially occupied molecular orbitals which, by virtue 
of Eq. (EH]), satisfy 



Vl k {t)JC{t ~ T)Vl m (T)(m\e s (r))dT 



Kk(t)n,n(t){m\es(t)) / £(t-r)dr. 

Jo 

Further, for times t 3> j^— the remaining time in- 
tegral is well approximated* by its asymptotic value, 
linit^oo J K,(t — r)dr = ft/ti oa d, and the equations of 
motion reduce to their Markovian limit: 



ih-{n\e s (t)) 



N 

E 

m— 1 



(n\Hg{t)\m) - i 



t 



coup -j 
^lcad 



where 



,{t)=Vl 1 {t)Vl m {t) 



N .rn 



He s (t)), 
(28) 



(t). (29) 



Equation (|28|) provides a realistic mechanism for elec- 
tron absorption and energy dissipation in the dynamics 
of a metal-molecule-metal system. The influence of the 
leads in the dynamics has been mapped into a negative 
imaginary (absorbing) potential. The projection opera- 
tors in the coupling, contained within r, vm , ensure that 
only those electrons with sufficient energy get absorbed, 
with proper conservation of the antisymmetry principle 
of the many-electron wavefunction. The equation is valid 
for times t ^> fi/ii ea d but can be applied for all times by 
turning on slowly the interaction between the molecule 
and the leads; a strategy that we adopt. 

Note that the field influences the dynamics directly by 
modifying the wire component, and indirectly by influ- 
encing the wire-lead coupling since T n ^ m {t) depends both 
on the nuclear dynamics and the radiation-matter inter- 
action. This characteristic introduced by the field has 
also been observed in an alternative description of the 
laser- induced dynamics of molecular wires [271 ] . Further 
note that this model neglects hole transport. Neverthe- 
less, since the SSH Hamiltonian is electron-hole symmet- 
ric, so that a subsequent inclusion of this contribution is 
just expected to double the resulting current. 



m— (n|es(*)) = Mo + a(u n+1 {t) - u n {t))\ (n + l|e s (t)> 

+ [-t + a(u„(t) - u„-i(t))} (n - l\e s (t)) 
+ \e\E(t) (na + u„(t))(n\e s (t)) 



ti 



''lead 



^2 T n>m (t)(n\es{t)}, 



(30) 

for n = 2, • • • , N — 1. At the boundaries of the oligomer, 
n = 1 and n = N, the orbitals satisfy Eq. ([50)) with- 
out the term that involves the amplitude of the orbital 
amplitudes outside the molecular domain. Equations ([6]) 
and (|30[) . together with the auxiliary quantity @, con- 
stitute a closed set of N(N + 2) coupled first-order differ- 
ential equations that are integrated using a Runge-Kutta 
method of order eight [4(|. The projective term in the 
dynamics T n , rn (t) is updated at every time step by di- 
agonalizing the system's electronic Hamiltonian H§(t) 
to obtain the instantaneous eigenstates |7s(0) defining 
the projection operator through Eq. ([T7]) . We project 
onto states that are above the Fermi energy, taken to 
be the zero reference energy. In other words, out of the 
N single-particle instantaneous eigenstates of Hg oc (t) we 
project onto the upper N/2 states. The metal-molecule 
coupling is turned on smoothly and slowly (in 10 fs) from 
preparation time. 

Electronic dephasing due to vibronic couplings is in- 
corporated by integrating the above equations of motion 
for an ensemble of initial conditions. The initial condi- 
tions are generated using the strategy detailed in Sec. 
IIC of Ref. Briefly, the starting optimal geometry of 
the oligomer is obtained by minimizing the total ground- 
state energy of the chain. Then, a normal mode analysis 
around this geometry is performed, yielding the ground 
state nuclear Wigner phase space distribution function in 
the harmonic approximation. By importance sampling 
this distribution, an ensemble of lattice initial conditions 
{u l (0) , p l (0)} is generated. The associated initial values 
for the orbitals {|e z (0))} are obtained by diagonalizing 
the electronic part of the Hamiltonian for each initial nu- 
clear geometry {u z (0)}. Each member i of the ensemble 
defines a quantum-classical trajectory and the set is used 
to obtain ensemble averages. 

Throughout this work we use the standard SSH param- 
eters for PA [l]| : a = 4.1 eV/A, K = 21 eV/A 2 , t = 2.5 
eV, M = 1349.14 eV fs 2 /A 2 and a = 1.22 A. In turn, we 
take the molecule and leads to be weakly coupled with 



coup 



lead 



= 0.1 eV. 



5. Final equations C. Dynamical observables 

Consider now the final equations of motion. The mi- The vibronic dynamics is characterized by the geomet- 

clear degrees of freedom satisfy trajectories determined ric and spectroscopic properties defined in Sec. IID of 
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Ref. l20j. In addition, an expression for the current that 
depends solely on molecular properties is required and is 
obtained in this section. 

The current entering into lead (3 =L, R is defined by 



J/3 



i ggg 

1 at 



(31) 



where pp is the number of electrons in lead (3. In terms 
of the reduced density matrix jp can be expressed as 



J/3 



= ^ |e| E (*(t)\lcl S Cn,s,H elec ]mt)h 



E 



_ (32) 

where we have used Eq. ([7])- For the Hamiltonian in 
Eq. (JSj> only the commutator with the lead-molecule cou- 
pling contributes to the current, 



IP = ftN E (*(*)l[4,.c n ,.»^s-/9]|*(*)>. (33) 

Introducing the projective coupling [Eq. (|22p] into the 
above expression, we have that 



2|e|t 



Jl = 



coup 



E MPl,m(*)P0,m(*)}, 



2|e|t. 



Jr. 



coup 



E M^,m(*)PJV+l,m(t)}. 



(34) 



m£S 



Equation (|34| indicates that the current entering into 
the leads depends on the spatial coherences between the 
lead boundary and every site in the molecular system. 
We desire, however, an expression that solely depends on 
molecular properties. For this we exploit Eq. ([T3"]l in its 
Markovian limit according to which, 



ricad m 



(35) 



so that 



Mead 
''lead 



(36) 



Substituting Eq. (|3"6"|) into Eq. (134(1 gives an expression 
for the current that is completely defined by molecular 
properties, 



2\e\ti 



^-lcad 

m.7i 



(37) 



as desired. Another useful quantity is the amount of 
charge that has been deposited in lead /? at a given time. 
This is defined by 



%{t) = / J>(r)d 



(38) 



Any rectification generated by an w + 2uj pulse results in 
qh - QR ^ 0. 



III. RESULTS AND DISCUSSION 

In Ref. important ways in which the electron- 
vibrational couplings influence the photoinduced dynam- 
ics of isolated PA chains were identified. The coupling 
introduces a significant exchange of energy between vi- 
brational and electronic degrees of freedom, broadening 
of the electronic transitions, pronounced changes in the 
mean single particle spectrum during and after the laser 
pulse, internal relaxation mechanisms and ultrafast deco- 
herence (in less than 10 fs). By contrast, in rigid systems 
no energy is transfered between the electronic and nu- 
clear degrees of freedom, and changes in the electronic 
spectrum are only due to Stark shifts induced by the 
laser field. 

Here we analyze the extent in which the vibronic cou- 
plings affect the ability of u> + 2u> lasers to induce cur- 
rents along molecular wires. Vibrational effects are made 
evident by contrasting the ensemble-averaged electron- 
vibrational evolution of the bridge with the dynamics of 
a single rigid trajectory that stays frozen at the optimal 
geometry. The wire is made rigid by arbitrarily increas- 
ing the mass of the (CH) groups by 10 6 . In this way the 
electron-ion interaction remains constant but the lattice 
moves a thousand times slower. Two different regimes 
of operation are considered: In Sec. IIII A[ the scenario is 
applied in its usual form where the field frequencies are 
tuned at or near an electronic resonance and electrons are 
promoted to the excited states through near-resonance 
photon absorption. In Sec. IIII Bl a novel mechanism is 
introduced in which the field is kept far off resonance and 
excitation is done through the dynamic Stark effect. As 
shown below, in the presence of vibrations the first sce- 
nario becomes inefficient while the Stark shift scenario 
remains robust. 



A. Currents through near-resonance excitation 

We investigate the dynamics and currents onset by the 
laser pulses specified in Table U for different laser frequen- 
cies and relative phases noted in the text. The param- 
eters chosen for the pulses are meant to encompass four 
illustrative cases: dynamics induced by weak and mod- 
erately strong pulses with time envelopes that are either 
comparable (10 fs) or long (300 fs) compared to the typ- 
ical electronic dephasing time. The laser frequencies are 
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TABLE I: Parameters and labels defining the femtosecond 
laser pulses used E(t) — exp( — (t — T c ) 2 /T^ v )(e UJ cos(wi+0 w ) + 
£2w cos(2tjt + <f)2u))- Here l2u> is the intensity of the 2u> com- 
ponent at maximum field strength. 



Label 


T c (fs) 


Tw (fs) 


€2^ (V A- 1 ) 


fu / £2w 


h u (W cm" 2 ) 


fl 


900 


300 


8.70 x 


2.82 


1.0 x 10 H 


f2 


900 


300 


4.00 x 10" 2 


2.82 


2.1 x 10 10 


f3 


50 


10 


8.70 x 10~ 3 


2.82 


1.0 x 10 9 


f4 


50 


10 


4.00 x 10 -2 


2.82 


2.1 x 10 10 



chosen to be at or near resonance so that, in this case, 
near-resonance photon absorption is the main source of 
electronic excitation. 

The system consists of neutral PA wires composed of 
20 sites and 20 tt electrons (L « 23 A) initially in the 
ground electron- vibrational state. An ensemble of 40,000 
initial configurations is propagated. This oligomer size 
was selected because it is representative of molecules that 
have been employed in constructing molecular nanojunc- 
tions The initial geometry of the chain consists of 
a perfect alternation of double and single bonds. The 
electronic structure is composed of 20 doubly occupied 
valence n orbitals and 20 empty n* states, separated by 
an energy gap of 2 A = 1.8 eV. A detailed account of the 
initial configuration, is provided in Ref. I20I . 




' 50 100 150 200 250 300 



t (fs) 

FIG. 1: Charge deposited in the left (solid lines) and right 
(broken lines) contact during and after photoexcitation with 
pulse f4. Here hui = 1.3 eV, §2^ — 2<j>u ~ (black lines, 
flexible wire; gray lines, rigid case). 

We begin by presenting two concrete examples of the 
currents induced on flexible and rigid PA wires by u> + 2uj 
pulses. Figures [1] and [2] show the charge deposited in the 
left and right lead during and after photoexcitation with 
a moderately strong 10 fs pulse (field f4 in Table [TJ) and 
a weak 300 fs pulse (field fl) resonantly coupling valence 
and conduction bands, respectively. Here, </>2 W — 4>u — 0. 
As can be seen, symmetry breaking is achieved while 
the system is being driven by the laser field. After the 
pulse, any remaining electrons with energy higher than 
the Fermi energy are absorbed symmetrically by the left 
and right contacts. At the particular frequency chosen 
we observe that the degree of symmetry breaking induced 




250 500 750 1000 1250 1500 1750 2000 



t (fs) 

FIG. 2: Total charge entering the left (solid lines) and right 
(broken lines) lead during and after photoexcitation with 
pulse fl with Hlo — 1.2 eV, 4>2ui — 2<^ = (gray lines, rigid 
wire; black line, flexible case). 
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FIG. 3: Frequency dependence of the net difference in charge 
deposited by pulse: (A) fl, (B) f2, (C) f3, (D) f4 with 2w - 
2(j>^ = 0. The black dots show the behavior of the flexible 
chain while the gray lines correspond to the rigid case. 



by the 10 fs pulse in the rigid and flexible chain is of the 
same order of magnitude. Note, however, that in the flex- 
ible wire considerably more electrons are excited across 
the energy gap than in the rigid example, due to the level 
broadening introduced by the wire's vibrations. Most of 
these electrons do not contribute to the rectification. As 
discussed below, this is a recurrent motif when applying 
the near-resonance lu + 2lu scenario to flexible systems. 
In the case of photoexcitation with a 300 fs pulse, and for 
the parameters chosen in Fig. [2] in rigid wires most of the 
photoexcited electrons participate in forming the dc cur- 
rent. In flexible wires, however, most excited electrons do 
not contribute to the net current. Although individual 
trajectories may exhibit large rectification effects, there 
is almost a perfect cancellation when the individual con- 
tributions are added together. The decoherence is simply 
too fast to allow maintenance of the rectification effect. 
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FIG. 4: Frequency dependence of the efficiency r\ of the rec- 
tification generated by pulse: (A) fl, (B) f2, (C) f3, (D) f4, 
with (j>2uj — 2(f>u> = 0. The black dots show the behavior of 
the flexible chain while the gray lines correspond to the rigid 
case. 
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FIG. 5: Net transport as a function of the relative phase 
4>2uj — 2</> CJ induced by pulse f4 with hui — 1.3 eV in the rigid 
(gray lines) and the flexible (open circles) chain. Panel (A) 
shows the net difference in charge deposited; (B) depicts the 
efficiency of the process. 



We now focus on the net asymptotic current induced 
by all of the laser pulses of Table Q] for different laser 
frequencies. The quantities of interest are the net dif- 
ference in charge deposited in the left and right leads 
— 9r, and the overall efficiency of the process r\ — 
(<Zl — ?r)/(?l + 9r)j where qp is the total charge de- 
posited at lead (3. Ideally one would like to find driving 
parameters that induce sizable currents in the nanojunc- 
tion with high efficiency. 

Figures [3] and [4] show <?l — 9r and rj after excita- 
tion with pulses fl-f4 for different laser frequencies, with 
4>2u) — 20 w = 0. Consider first the rigid chain results. For 
rigid chains the degree of control and the net induced cur- 
rents are considerable. The effect is enhanced by tuning 
the frequency components of the laser near one of the 
system's electronic transition frequencies. At selected 
frequencies the process can achieve efficiencies as high 
as ~ 60%. Note, however, that for exact on-resonance 
driving frequencies the efficiencies are low due to com- 
peting multi-photon processes that populate the excited 
states but do not contribute to the currents. By applying 



stronger pulses (Fig. [3)3) it is possible to exploit higher 
order multiphoton processes to induce rectification, re- 
sulting in a complicated frequency dependence to the 
control map. When short pulses are employed (Fig. 
D) the sharp features observed in the long-pulse control 
are no longer resolved. For rigid wires, higher efficiencies 
are achieved by using long weak pulses at frequencies that 
minimize satellite photoexcitations. When employing ei- 
ther short or very strong pulses, satellite absorptions, 
arising from either competing multiphoton processes or 
parasite absorption at frequencies contained within the 
bandwidth of the pulse, are unavoidable. 

The inclusion of vibrations in the dynamics substan- 
tially modifies the physical picture. As shown in Figs. [3] 
and 21 lattice fluctuations smear out the resonances ob- 
served in rigid wires resulting in control maps with broad 
features only. Although the scenario is successful in gen- 
erating phase controllable currents (see, e.g., Fig. [5]) with 
magnitudes that are often larger than the ones observed 
in rigid wires, the efficiency of the process is invariably 
very low irrespective of the driving frequency, laser inten- 
sity or bandwidth (Fig. [4]). That is, most photoexcited 
electrons do not participate of the currents and the en- 
hancement of the effect is due to increased light absorp- 
tion of the flexible wire relative to the rigid case. The 
highest observed efficiencies, obtained for moderately 
strong short pulses (Figs. GMP), are merely r\ 0.04. 
These results clarify the previously unexplained origin 
for the enhancement of currents b y in creasing vibronic 
coupling strength observed in Ref. [301 ] . 



B. Currents through the dynamic Stark effect 

Thus far our approach has been to excite electrons 
across the energy gap of the wire by tuning the frequen- 
cies of the field near a resonance of the system. We have 
seen, however, that this technique is fragile to decoher- 
ence processes since it relies on creating coherent super- 
position states. Lattice fluctuations make the rectifica- 
tion extremely inefficient and, in the best of cases, only 
about 4% of the electrons that are photoexcited partic- 
ipate in the net current. Hence, in order to induce siz- 
able currents considerable energy from the field must be 
dumped into the nanojunction, compromising its struc- 
tural integrity. Further, employing faster or stronger 
pulses is not very helpful since they introduce other un- 
desired satellite channels in the control and do not ap- 
preciably overcome the deleterious effects introduced by 
the vibrations. 

In this section we introduce an alternative mechanism 
that is remarkably robust to electron-vibrational cou- 
plings, survives in the presence of decoherence and ther- 
mal effects, and is able to induce large currents in molec- 
ular wires with efficiencies > 90%. Instead of promoting 
electrons to the conduction band through near-resonance 
photon absorption, we work far off resonance and em- 
ploy Stark shifts to nonadiabatically couple the ground 
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and excited electronic states. Phase-controllable symme- 
try breaking is achieved by exploiting the difference in 
the intensity of u> + 2lo fields for positive and negative 
amplitudes. 

Although this mechanism is applicable in medium sized 
chains (N w 20) it is best suited for long oligomers. 
Longer wires have a smaller energy gap, facilitating the 
coupling of ground and excited states through Stark 
shifts before photon absorption becomes dominant. For 
this reason we consider perfectly dimerized neutral PA 
oligomers with N = 100 sites (L w 120.8 A) originally in 
the ground state configuration, and propagate an ensem- 
ble of 1000 initial configurations. The electronic struc- 
ture consists of N/2 doubly occupied tt (valence) orbitals 
and A^/2 empty tt* (conduction band) states, separated 
by an energy gap of 2 A = 1.3 eV which is significantly 
smaller than the 1.8 eV observed by 20 site oligomers. 
The average total energy of the ensemble is ~ 4.4 eV 
higher than its rigid counterpart due to zero-point fluc- 
tuations of the lattice in the ground electronic surface. 



1. Complete control 

We now illustrate the essentials of the phenomenon 
and show how under realistic conditions one can achieve 
almost complete control over the currents even in the 
presence of significant electron-phonon couplings. Con- 
sider the clcctron-vibrational dynamics of the wire under 
the influence of an u> + 2ui field of the form 



E(t) 



n=l 

exp ( 
1 



u cos(nLdt + nw )x 
-(t-t on ) 2 /T*) for t<t 



(39) 



for t on < t < toS 



exp (-(< - t oS y/T£) for t > t oS 



The field is smoothly turned on and off in T w — 100 fs 
(t on = 300 fs; t ff = 700 fs) and has constant amplitude 
for 400 fs. The frequency hiu = 0.2A = 0.13 eV is chosen 
far off-resonance from the system's interband transition 
frequencies so that Stark shifts, and not photon absorp- 
tion, dominate the dynamics. The field amplitude used 
is = 6.1 X 10~ 3 V A -1 with e w = 2e 2w , which corre- 
sponds to an intensity I 2uJ ~5x 10 8 W cm~ 2 . 

Figure [5] shows the field, the currents, and the light- 
dressed electronic structure of the flexible chain, aver- 
aged over all trajectories, when the relative phase of the 
pulse is <t>2u> — 2</>o; = 0. Due to the highly polarizable 
nature of 7r-conjugated systems, the single-particle spec- 
trum displays considerable Stark shifts, even at moder- 
ate field intensities; the effect being stronger for states 
near the edges of the valence and conduction band. In 
this way, the dynamic Stark effect effectively reduces the 
energy gap of the oligomer, causing frequent crossings 
between the ground and excited electronic states in indi- 
vidual trajectories. Since the wire's ground state is non- 
degenerate and of definite parity, the lowest-order contri- 
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FIG. 6: Time dependence of (A) the electric field, (B) the 
current entering the left and right leads, and (C) the instanta- 
neous light-dressed orbital energies for states near the energy 
gap for a 100-site flexible wire under the influence of the field 
in Eq. (|39[l with <p2uj — 2<j>u = 0. Note the bursts of charge 
deposited in the right lead when the field bridges the energy 
gap. Here j'l is so small that is barely visible. 




FIG. 7: Time dependence of (A) the electric field; (B) the 
current entering the left and right leads and; (C) the instanta- 
neous light-dressed orbital energies for states near the energy 
gap for a 100-site flexible wire under the influence of the field 
in Eq. (|39[) . with 02a; — 20^ = n/2. Note that in this case 
the bursts of charge deposited alternate between the left and 
right contact. 



bution to the Stark effect is quadratic in the field. Hence, 
when |£7(t)| is maximum the energy gap acquires acquires 
its minimum value. At the crossing times population is 
transferred from the valence to the conduction band and 
bursts of charge are deposited in the leads. 

Note that for — 20^ = almost all excited elec- 
trons are deposited in the right contact only. Symmetry 
breaking arises due to the difference in the maximum 
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fact, the currents observed in the flexible case are al- 
ways larger than in the rigid example. However, in rigid 
wires the mechanism can exhibit perfect efficiencies. The 
mechanism is also expected to be resilient to thermal ef- 
fects. Increasing the temperature will introduce addi- 
tional broadening of the electronic levels and, as pointed 
above, the Stark rectification is robust to this class of 
effects. 



FIG. 8: Phase dependence of (A) the net rectification, and (B) 
efficiency of the process when the field in Eq. (|39|l is applied to 
flexible and rigid 100-site PA wires (black dots, flexible wire; 
gray dots, rigid case). Note the assistance of phonons in the 
rectification. 



|.E(i)| for positive and negative amplitudes exhibited by 
the field. Even when E(t) has a zero temporal mean, 
it consists of narrow peaks with large \E(t)\ for positive 
amplitudes, and shallow and broad features when E(t) 
is negative (Fig. [6^). For this reason, the Stark effect is 
only sufficiently strong to close the energy gap when the 
field has a positive amplitude. Thus, transfer of popula- 
tion to the conduction band and absorption of electrons 
by the leads always occur when the laser is pointing at a 
particular and the same direction, in this way inducing 
directed transport in the system. 

The phenomenon depends intimately on the relative 
phase. For instance, for the case of 4>2u — 20 w = k/2 
shown in Fig. [7] the electric field exhibits equal inten- 
sity for positive and negative amplitudes. Since the field 
changes sign between consecutive interband couplings, 
the bursts of charge deposited alternate between the left 
and right contact, and no net current is induced. 

Figure [8] shows the net difference in charge deposited 
in the left and right leads after the pulse is over for differ- 
ent laser phases, as well as the efficiency of the process. 
For comparison purposes the plot also includes results 
obtained in an equivalent but rigid system. We first note 
that the mechanism is robust to decoherence effects due 
to coupling to the vibrational degrees of freedom as well 
as satellite contributions due to parasite multiphoton ab- 
sorption. In fact, 90% of the excited electrons can par- 
ticipate in the net current. This should be contrasted 
to the extremely low efficiencies achieved before through 
the multi-photon absorption mechanism. Further, the 
sign and magnitude of the effect can be manipulated by 
varying the laser phases. By changing the relative phase 
by 7r the magnitude of the effect stays the same but the 
direction of the rectification is reversed. 

In the flexible wire the rectification exhibits an almost 
sinusoidal dependence on <f>2uj — 20 w . By contrast, in 
rigid wires for certain range of phases no currents are 
induced since the maximum field amplitude is not large 
enough to couple valence and conduction bands. Hence, 
in this range the currents observed in the flexible wire 
are phonon-assisted; i.e. the level broadening introduced 
by the vibrations permit the nonadiabatic coupling. In 



2. Possible satellite contributions 

The simulations presented above exemplify how almost 
complete control over the electronic dynamics in the pres- 
ence of vibronic couplings can be exerted in molecular 
wires. However, it is interesting to point out possible 
sources of satellite contributions that reduce the effi- 
ciency of this rectification mechanism. For definitiveness 
consider the wire's dynamics under the influence of an 
lu + 2lu Gaussian pulse of the form 

2 

E{t) = e -^- T ^l T « tnu cos(nut + (40) 

centered at T c = 900 fs, with temporal width Tw = 300 
fs, and amplitudes t2u> = 8.7 x 10~ 3 V A -1 with e w = 
2e2u, which corresponds to an intensity of l2 U = 10 9 
W cm -2 , and driving frequency hui — 0.46A = 0.3 eV. 
With respect to the field in Eq. (|3"9"|) . this laser profile 
has a higher maximum intensity and blue shifted central 
frequencies. Figure [9] shows the currents onset by this 
laser pulse when the relative phase is 4>2u> — = 0, and 
Fig. [TO] the total charge deposited during the dynamics. 
The field generates net currents along the system, with 
the right contact accepting 1.25 more electrons than the 
left one. However, in this case only 9% of the excited 
electrons participate in the net current. 

There are two sources of satellite contributions that 
reduce the efficiency of the scenario. The first is that 
some interband coupling occurs when the amplitude of 
the field is negative, contrary to the case described in 
the previous section. That is, the shallow negative am- 
plitude of the field is sufficient to close the energy gap in 
some trajectories and deposit electrons in the left lead. 
Second, since the field frequency is larger than the one in 
Eq. (|39[) there is considerable excitation originating from 
photon absorption and not from Stark shifts. This leads 
to a background signal in jp (t) that does not contribute 
appreciably to the rectification and reduces the efficiency 
of the process. It is not difficult to suppress these two 
satellite effects by changing field parameters. Parasite 
multiphoton absorption can be avoided by further detun- 
ing the frequency of the field from the interband transi- 
tion frequencies. Undesired closing of the energy gap 
can be circumvented by employing a laser profile with 
controlled and reduced intensity. Highly efficient current 
generation, such as that seen in the previous example, 
results. 
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FIG. 9: Time dependence of (A) the electric field, and (B) the current entering the left and right leads for a 100-site flexible 
wire under the influence of the field in Eq. (|40[) with 4>2^ — 2<^„ = 0. The insets detail the field and currents from t — 900 fs to 
1000 fs. 
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FIG. 10: Total charge deposited in the left §l (black line) and 
right qn (gray line) contact during the dynamics of a 100-site 
flexible wire under the influence of the field in Eq. (|40|l . with 
4>2^ — 2<f>u = 0. The dotted lines provide an estimate of the 
sampling errors. 



IV. CONCLUSIONS 

Here we have studied the possibility of inducing ultra- 
fast currents using lu + 2lo fields along molecular wires 
with significant electron-vibrational couplings. Two 
possible mechanism for laser-inducing rectification were 
identified. In the first one, electrons are excited across 
the energy gap of the wire by tuning the frequencies of 
the field near an electronic resonance of the system. By 
applying lasers with both lu and 2cu components asym- 



metry is generated in the momentum distribution of the 
photoexcited electrons. In the second mechanism, the 
lu + 2lu field is kept far off-resonance and population is 
transferred to the excited states by means of Stark shifts 
that bridge the energy gap of the oligomer and, in this 
way, couple ground and excited electronic states. Sym- 
metry breaking is achieved by exploiting the difference 
in the field's intensity for positive and negative ampli- 
tudes. This permits transferring population across the 
energy gap only when the field is pointing at a particular 
direction, thereby breaking the spatial symmetry of the 
system. In both cases the magnitude and direction of the 
net currents is controllable by varying the relative phase 
of the two frequency components. 

The vibronic couplings in the wire have very differ- 
ent effects in these two cases. In the first one, ultrafast 
decoherence processes induced by the wire's vibrations 
make the rectification effect extremely inefficient. While 
in rigid wires the maximum efficiency observed is ~ 60%, 
in flexible chains less than 4% of the photoexcited elec- 
trons participate in the current. This trend holds irre- 
spective of the driving frequency, field strength or width. 
At selected frequencies the vibrations can enhance the 
magnitude of the effect, but not the efficiency, due to in- 
creased light absorption through level broadening. This 
observation sheds light on thepreviously unexplained ob- 
servations reported in Ref. [301 ] where the magnitude of 
the photoinduced currents was shown to increase with 
the strenght of the vibronic couplings. Note, however, 
that in order to generate sizable currents in this regime 
a considerable number of electrons must be excited, un- 
dermining the structural stability of the system. 

In contrast to the above mechanism, the off-resonance 
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Stark mechanism introduced here is impervious to vi- 
bronic couplings. Even when the lattice dynamics may 
introduce parasite multiphoton absorption or undesired 
closing of the energy gap, these problems are easily over- 
come by red shifting the driving frequency and using 
pulses of lowered and controlled intensity. Simulations 
show that the mechanism can actually induce currents 
with efficiencies as high as 90%, in spite of the significant 
electron-phonon couplings. 

We note that the actual implementation of these 
schemes poses substantial experimental challenges. The 
laser beam needs to be focused on the wire and the setup 



must control possible side effects that may arise from 
the illumination, such as heating, or processes involv- 
ing excitation of the metal surface. Recent experimental 
developments [U [42|, HH suggest that this may be pos- 
sible. We hope that emerging applications, like the one 
presented here, stimulate further experimental progress. 

Last, from the simulations presented herein and a re- 
cent experiment |44| . it is clear that the dynamic Stark 
effect is a promising tool to exert laser control that de- 
serves further exploration. The success of this approach 
relies on the fact that it does not exploit the fragile prop- 
erties of superposition states. 
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